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We discuss the low-energy limit of three-orbital Kondo-lattice and Hubbard models describing t2 g 
orbitals on a triangular lattice near half filling. We analyze how very flat bands with non-trivial 
topological character, Chern number C = 1, arise both in the limit of infinite on-site interactions 
as well as in more realistic regimes. Exact diagonalization is then used to investigate fractional 
filling of an effective one-band spinless-fermion model including nearest-neighbor interaction V; it 
reveals signatures of fractional Chern insulators (FCIs) for several filling fractions. In addition 
to indications based on energies, e.g. flux insertion and fractional statistics of quasiholes, Chern 
numbers are obtained. It is shown that FCIs are robust against disorder in the underlying magnetic 
texture that defines the topological character of the band. We also investigate competition between 
FCI states and a charge density wave (CDW) and discuss particle-hole asymmetry as well as Fermi- 
surface nesting. FCI states turn out to be rather robust and do not require very flat bands, but 
can also arise when filling or an absence of Fermi-surface nesting disfavor the competing CDW. 
Nevertheless, very flat bands allow FCI states to be induced by weaker interactions than those 
needed for more dispersive bands. 



I. INTRODUCTION 

Quantum liquids are among the most sought after 
states of matter. One celebrated class of such quantum 
liquids, the fractional quantum-Hall (FQH) states, 1-3 
has been a locus of scientific attention for almost three 
decades. Initially introduced to explain the fractional 
quantum-Hall effect (FQHE) observed in semiconductor 
devices, 4 they describe interacting electrons constrained 
in two dimensions and subject to a strong perpendicu- 
lar magnetic field. Fractionally filling a magnetic Lan- 
dau level (LL) then yields an incompressible FQH liquid 
and gives rise to precisely quantized Hall conductivity. 
In addition, the quasiparticles of these states obey any- 
onic statistics, 5,6 which can be either Abelian or non- 
Abelian, the latter fulfilling an essential condition for 
fault-tolerant quantum computation. 7 

An alternative route to quantum Hall liquids is via 
tightly bound electrons moving in a magnetic texture. 
The simplest such magnetic texture is a uniform mag- 
netic field perpendicular to the plane of the system. For 
non-interacting, tightly bound electrons, this gives rise to 
an integer quantum-Hall (IQH) state and leads to Hofs- 
tadtcr's fractal energy spectrum, 8 whereas in a fraction- 
ally filled tight-binding model of electrons interacting via 
a screened Coulomb repulsion, the realization of a FQH 
state is possible. 9 Haldane noted 10 that a magnetic field 
is not the only viable path to FQH state and introduced a 
honeycomb lattice model with complex hoppings, which 
break time-reversal invariance, and thus defined the first 
integer Chern insulator (CI), where the total magnetic 
field through the system averages to zero. In the mean- 
time, mechanisms leading to complex hoppings with the 
necessary properties have been identified. One realiza- 
tion may be found in strongly spin-orbit coupled semicon- 
ductor materials that are ferromagnetically ordered. 11,12 
Another possibility arises through the coupling of itin- 
erant electrons to localized magnetic moments 13 , for ex- 



ample the Kondo-latticc model on the triangular lattice 
supports a non-trivial magnetic texture, which induces 
an integer-quantized Hall conductivity of the itinerant 
electrons. 14 

A considerable body of recent research has addressed 
the question whether the lattice counterpart of the FQHE 
can be observed when Chern bands are fractionally filled 
and electrons are interacting. Several numerical studies 
using exact diagonalization techniques have convincingly 
established the existence of Laughlin-series states in an 
number of different models on various lattices. 15 " 21 These 
systems have since been called fractional Chern insula- 
tors (FCI). Very recent work has reported FCI states be- 
yond Laughlin fractions. 22-24 In addition to reproducing 
the known FQHE on a lattice and potentially at higher 
temperatures, CIs also offer the intriguing possibility of 
a Chern number higher than one, a departure from the 
analogy with LLs. Recent studies have explored this di- 
rection by constructing models that have higher Chern 
numbers and studying possible FCI states. 25-29 From the 
analytical side the problem of FCI states has been ap- 
proached by a careful study of emergent translational 
symmetries 30 and many-body trial wave functions. 31-33 
Others have examined the algebraic properties of the den- 
sity operators projected onto one Chern band and made 
a comparison with the Girvin-MacDonald-Platzman al- 
gebra that is satisfied by the lowest LL density operators 
in the continuum FQHE. 30,34-37 

In such studies of FCIs, one typically starts by adding 
inter-site interactions to topologically nontrivial but non- 
interacting CI models. As potential realizations of the 
non-interacting and nearly flat "parent" bands, cold 
atoms, 16 oxide heterostructures, 38 strained graphene, 39 
and strongly correlated multi-orbital models for layered 
oxides 22 ' 40 have been proposed. In the present paper, 
we build upon this last approach and thus focus our at- 
tention on a strongly correlated model on the triangular 
lattice. It was shown that on the mean-field level and 
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near half filling, a magnetically ordered CI with a very 
flat band emerges. Doping this nearly flat band to frac- 
tional fillings was shown to give rise to FCI states within 
the framework of an effective single-orbital model. Here, 
we give a detailed account of the mapping onto the effec- 
tive single-orbital model and show how such topologically 
nontrivial and nearly flat bands emerge in Kondo-lattice 
and Hubbard models both for the limit of infinite onsitc 
interactions and for more realistic intermediate interac- 
tion strength. 

We then provide extensive numerical evidence for the 
existence of FCI states in this model by using exact di- 
agonalization. We discuss the robustness of FCI states 
against disorder originating from single-site defects in the 
magnetic ordering, which will always occur in realistic 
situations. Another very relevant issue in the context of 
FCI states is their competition with other phases, 24 ' 29 
e.g. symmetry-broken states such as a charge-density 
wave (CDW). We study this issue by using a method that 
does not project onto the nearly flat Chern band and thus 
keeps the effects of dispersion. The model considered in 
this work allows for a careful study of the competition 
between finite dispersion and interactions. We map out 
a phase diagram for filling fractions v — 1/3 and v = 2/3. 
In the latter case, filling permits a commensurate CDW, 
whereas in the former it does not. The CDW is favored 
by Fermi-surface (FS) nesting, and we accordingly find 
the FCI at v = 2/3 to be far more stable when the bands 
are poorly nested. 

This paper is organized as follows: in Sec. II, we 
extend the discussion of flat and topologically nontriv- 
ial bands arising in t-2 g models and derive the effec- 
tive model. Numerical results on the model are pre- 
sented in Sec. Ill, where we first focus on information 
obtained from cigenenergies, namely gaps and flux in- 
sertion (Sec. Ill A) and fractional statistics of charged 
excitations in FCI states based on a recently introduced 
state-counting argument (Sec. Ill B). 30 We then add in- 
formation obtained from the eigenstates: the Chern num- 
ber, which allows us to address the impact of impurities 
in Sec. Ill D and the static charge-structure factor, which 
allows us to discuss the competition with the CDW at fill- 
ing v = 2/3 in Sec. HIE. We conclude with some remarks 
summarizing the main points in Sec. IV 



II. TOPOLOGICALLY NONTRIVIAL AND 

NEARLY FLAT BANDS IN STRONGLY 
CORRELATED MULTI-ORBITAL SYSTEMS 

In Ref. 40, it was shown that both e g and t2 g or- 
bital manifolds in octahedral coordination can reduce the 
bandwidth of topologically nontrivial bands. In partic- 
ular, this was discussed for the spin-chiral phase aris- 
ing in Kondo-lattice models 14, 41-43 at quarter and three- 
quarter fillings. While the flat band of interest mixes 
both the 3z 2 — r 2 and the x 2 — y 2 orbitals in the e g case, it 
is dominated by a particular orbital state in the ti g man- 
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FIG. 1. (Color online) Illustration of chiral spin pattern and 
effective model including the magnetic texture as phase fac- 
tors in hoppings. 



ifold, the a\ g state, which allowed the straightforward 
mapping onto an effective one-band model. 22 In this sec- 
tion, we discuss how the nearly flat a\ g band arises due to 
effective longer-range hoppings. The section is divided in 
subsections, where we first discuss the simplest scenario 
of infinite Hund's rule coupling in Sec. II A, and later the 
relation to the more realistic multi-orbital model with 
finite interactions in Sec. II B 



A. Effective bands for a Kondo-lattice model with 
infinite Hund's rule coupling 

The impact of the orbital symmetries on the hopping 
is contained in orbital- and direction- dependent matrix 
elements & , where a. and (3 denote orbitals (xz, yz, 
and xy) and the direction. We only take into account 
nearest-neighbor hopping elements, which are for i 2 g or- 
bitals on the triangular lattice given by 
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for NN bonds along the three directions ai, a 2 , a3. 
Direct orbital overlaps give orbital-conserving elements 
t ddl while oxygen- mediated processes lead to interorbital 
to- 44 ' 45 If the width of a triangular layer made of octa- 
hedra is compressed (extended), the energy of the highly 
symmetric orbital state \a\ g ) = {\xz) + \yz) + \xy})/\/3 
is raised (lowered) with respect to the remaining orbital 
doublet (e' g ). This energy shift can be written as 



H JT = -A JT (n e + 
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and depends on the Jahn- Teller effect as well as the lat- 
tice. 45 

The interplay of these orbital symmetries with the 
Berry phases arising from a spin-chiral phase was stud- 
ied in the limit of infinite Hund's rule coupling to classi- 
cal localized spins, corresponding to the double exchange 
model. In this case, electrons become spinless fermions, 
as only the spin projection parallel to the localized spin is 
relevant, it is designated as 'up'. Hopping between sites 
i and j = i + a.k (connected by a bond along lattice di- 
rection afc) is modulated by the Berry's phase mJJ 46 and 
becomes 
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where the classical spin at site i is parametrized by angles 
9i and fa , with the shorthand notation Cj = cos 6i and 
Sj = sin^. For the chiral spin pattern in Fig. 1, the 
Berry phases between the four sites of the magnetic unit 
cell can be parametrized as 
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It turns out that these hoppings can in fact be written 
in a two-site unit cell (containing sites 1 and 2) due to 
an internal symmetry of the four-site pattern. 14 . Using 
this two-site electronic unit cell as a basis, phases utt 
between NN unit cells are proportional to the three Pauli 
matrices along the three lattice direction, with the two 
states referring to the site within the unit cell, not the 
electron spin. 

The non-interacting spinless three-orbital model with 
hoppings given by Eqs. (1), (2) and (4) can be easily 
solved in momentum space and one finds that large |Ajx| 
strongly reduces the dispersion of one subband. 40 This 
can also be seen in Fig. 2(a), which shows the one-particle 
energies obtained on a cylinder. Figure 2(a) also reveals 
the edge states crossing some gaps, indicating the topo- 
logically nontrivial nature of these bands. Calculating 
Chern numbers C corroborates this and gives C = ±1. 
The band flatness can be expressed in terms of a figure 
of merit 
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where Ajt and A c are the two gaps separating the nar- 
row band of interest from the other orbitals and from the 
subband with opposite Chern number and W is the width 
of the narrow band. As can be pointed out, 40 the lower 
subband can here become very flat, and as can be seen in 
Fig. 2(b), the flatness can be further improved by going 
to larger crystal fields and can reach values M w 28. 

As these very flat bands can be achieved for large sep- 
aration Ajt between the a\ g and e' g states and as the 
band of interest then has almost purely a\ g character, it 
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FIG. 2. (Color online) Flat lower chiral subband in 
the Kondo-lattice model with infinite Hund's rule coupling 
(double-exchange model), (a) Shows the one-particle ener- 
gies of three t? g orbitals coupled to localized spins, where the 
latter form a spin-chiral phase on a triangular lattice, 14,41-43 
see Fig. 1. The system is a cylinder, i.e., periodic bound- 
ary conditions along y-direction and open boundaries along 
x. The horizontal axis is the momentum in the direction with 
periodic boundaries. The gaps Ajt and A c denote the gaps 
due to crystal-field splitting and to the chiral spin state, (b) 
shows the figure of merit M, see Eq. (5), for the lower ai g sub- 
band. The curves for crystal-field splittings Ajt = 4, 4.5, and 
5 were already given in Fig. 3(b) of Ref. 22 and are repeated 
here for convenience. 



is natural to assume that one should be able to capture 
the most relevant processes with an effective a\ g model. 
This is in contrast to the situation starting from e g or- 
bitals, where one finds intermediate Ajt to be optimal. 40 
In that case, the nearly flat bands can only be obtained 
if both orbitals contribute weight and one cannot easily 
reduce the situation to a one-band system. 

We now go over into the (ai g , e' g lr e' g 2 ) basis, which is 
done via 45 
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The transformed hopping matrices T ai are then obtained 
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from Eq. (1) as 
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where St = tdd — to and u> = e l27r / 3 . The direct hopping 
of an a lg electron into the a lg state of a NN site is thus 
isotropic (as expected for a\ g symmetry) and given by 



ti,j = h = ( 2t o + tdd)/3 u 
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However, the spin-chiral phase with only NN hoppings 
and infinite Hund's rule would not have particularly flat 
bands, but only reaches M w 2. 

The impact of the e' g levels can then be taken into ac- 
count in second-order perturbation theory. This includes 
processes where a hole hops from the a\ g orbital at site 
i to an e' g state at j and back again to an a lg state at 
a third site i' , which may or may not be the same as i. 
The denominator of these terms is the crystal-filed en- 
ergy Ajt and the numerator is obtained from the prod- 
ucts f? b f!> a +f? c f!f a (with a designating a lg and 6, c the 
e' states). In order to evaluate the second-order hopping 



FIG. 3. (Color online) Figure of merit M, see Eq. (5), for finite 
Hund's-rule coupling JKondo/io and Ajt = 6f . The bands 
designated as "upper" and "lower" refer to the two subbands 
of the dig states with spin parallel to the localized spin, which 
are separated by the gap opening in the spin-chiral phase, see 
Fig. 2(a). 



orbital double-exchange model, i.e., for infinite Hund's 
rule coupling to some additional localized spins, this sec- 
tion will discuss finite Hund's rule. Figure 3 shows the 
figure of merit for the band flatness Eq. (5) for a few sets 
of hopping parameters and for Ajt = 6*o depending on 
Hund's rule coupling to the localized spin. Hund's rule 
coupling is given by 



between sites i and i' ', these orbital hoppings have to be i?Kondo = — ^Kondo SjSi Q = — JKondo ^^(^q 



multiplied by the product of the Berry phases uj^ and 
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W- \, from Eq. (4) for all paths connecting i and i' via one 
intermediate site j ^= 



Due to destructive interfer- 
ence, processes connecting NN and next-nearest neighbor 
(NNN) sites cancel while effective third-neighbor hop- 
ping, where there is only one path, remains. Since third- 
neighbor spins in the chiral phase are always parallel, 
the total Berry phase of this process is 1 in all direc- 
tions, however, the hopping via a spin of different orien- 
tation in the middle reduces the hopping amplitude by 



D/2 
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where Si is the electron spin on site i and in orbital a, Si 
is the localized spin, and n\ (n^) is the electron density 
at site i in orbital a and (anti-) parallel to the localized 
spin. The latter is assumed to be large and thus treated 
as a classical spin of length 1. As electrons can now 
also be antiparallel to the core spin, even though this is 
energetically disfavored, additional Berry phases 
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A third-neighbor hopping oc '^2 i cos 2ka^ turns out to 
have almost the same dispersion as the chiral subbands 
and can consequently almost cancel it in one subband. 
As its strength can be tuned by tdd and Ajt, very flat 
subbands can be achieved, see Fig. 2(b). 



Impact of the upper Kondo/Hubbard band: 
Flat bands for finite interactions 



While the previous section has illustrated how one can 
understand the occurrence of nearly flat bands in a three- 



can come into play, with a = ±1 and a = —a denoting 
electron spin parallel and anti-parallel. As can be seen 
in Fig. 3, the upper subband of the a\ g sector can now 
become nearly flat. (For JKondo |Ajt|, one can of 
course still find flat lower subbands, as discussed above.) 

The flatness of the upper subband can be explained by 
similar effective longer-range hoppings in second-order 
perturbation theory, this time also taking into account 
intermediate states with an electron in the upper Kondo 
band, i.e. with antiparallel spin. These additional terms 
can go either via the a\ g or via the e' g orbitals and involve 
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has to sum over all possible intermediate sites j and finds 
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The NN, NNN and third-neighbor hoppings are here de- 
noted by ii, and t 3 . E x — A JT , E 2 = JKondo, and 
i?3 = JKondo — Ajt give the excitation energies of the 
intermediate states with (i) a hole in the e' g states with 
spin parallel, (ii) an electron in the a\ g states with spin 
anti-parallel and (iii) an electron in an e' g state with spin 
anti-parallel. Like the bare NN hopping, these effec- 
tive hopping-matrix elements acquire an additional Berry 
phase uJJ, in the Hamiltonian, which only depends on the 
relative orientation of spins on the initial and final sites. 
NNN hopping t 2 via the upper Kondo band does not drop 
out, and NN hopping becomes somewhat renormalized. 

The flat chiral subbands that have been observed in 
a three-orbital t 2g Hubbard model on the triangular 
lattice 22 arise in situations similar to the finite- JKondo 
scenario, and a mapping to the Kondo-lattice model 
was shortly mentioned in the supplemental material of 
Rcf. 22. It amounts to noting that the e' g levels are 
half-filled and far from the Fermi level and thus form 
a localized spin, while the a\ g states are closer to the 
Fermi energy and provide the itinerant electrons of the 
Kondo-lattice model. For infinite onsite interactions U 
(Coulomb interaction) and J (Hund's rule coupling), 
the three-orbital Hubbard model maps to the double- 
exchange model with infinite JKondo for large |Ajt| 
model and fillings close to half filling. 

However, the ground state is then a ferromagnct, as 
the Kondo-lattice model requires either finite JKondo 41 
or additional AFM inter-site superexchange 42 to support 
a spin-chiral state. At finite onsite interactions, second- 
order processes similar to the effective hoppings discussed 
here arise that involve e' g electrons and do not mediate 
hopping but instead generate magnetic superexchange 
interactions. This mechanism effectively provides the 
needed AFM intersite superexchange and the spin-chiral 
becomes stable for wide parameter regimes, including 
ranges supporting nearly flat upper chiral subbands. 22 



C. Effective Model 

As we have discussed in the previous Sec. II B, the 
most realistic route to nearly flat bands with nontrivial 
topology on the triangular lattice arises at finite Hub- 
bard/Kondo coupling, where effective second-neighbor 
hopping is generated in addition to NN and third- 
neighbor terms. However, we have found that all essential 
features of the band structure can be captured by using 
just NN and third-neighbor hoppings (see Figure in the 




FIG. 4. (Color online) (a) Energy dispersion ek for t' /t — 0.2. 
The inset shows the path taken through the first Brillouin 
zone, (b) Flatness ratio M, see Eq. (5), as a function of t' /t. 
The flatness ratio has been calculated from the dispersion 
along the high-symmetry directions shown in the inset on the 
left. 



supplemental material of 22). In the interest of simplic- 
ity, we consequently adopt here the kinetic energy 

#kin(k) = 2t ^2 a3 cos + 2t> a ° cos 2ka J ' ( 15 ) 

where &j (j = 1,2,3) denote the unit vectors on the trian- 
gular lattice, and t and t' are the NN and third-neighbor 
hopping, which can be related to Eqs. (8) and (9) for the 
double-exchange scenario, to Eqs. (12) and (14) for finite 
onsite interactions, or which can be taken as fit param- 
eters. Pauli matrices er- 7 and unit matrix er° refer to the 
two sites of the electronic unit cell in the chiral state. 14 
The unit cell and the topologically non-trivial bands are 
due to the symmetry breaking involved in the underlying 
magnetic order. The dispersion of Eq. (15) is 
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From now on, we use the effective NN hopping t as unit 
of energy; the band flatness can then be tuned by varying 
the ratio t' /t. 

The longer range hopping t' determines the flatness of 
the bands of /Jkin, which can be expressed by the figure of 
merit M, see Eq. (5). Figure 4 shows M depending on t', 
and one sees that ratios M > 20 can be reached for t' « 
0.2. Such flatness ratios can reasonably be achieved in 
the low-energy bands of a strongly correlated t 2g system 
on a triangular lattice. 22 Changing the sign of t' simply 
mirrors the dispersion vertically, i.e. it is then the upper 
band that becomes nearly flat. When going away from 
maximal M, the bands for smaller and larger t' differ 
qualitatively: For t' < 0.2, the Fermi surface (FS) at 
some fillings is almost perfectly nested. We are going to 
discuss the impact of these differences in Sec. HIE. 

NN Coulomb interaction 



(17) 
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is added to the kinetic term Eq. (15), giving the total 
Hamiltonian H = ifkm + flint- While the spin-chiral 
phase providing the nontrivial topology can only be ex- 
pected to remain stable for not-too-large V and doping 
of the flat band, we are going to study a variety of filling 
and interaction ranges to obtain a comprehensive picture 
of the model on finite-size systems. 



III. RESULTS 

The interacting Hamiltonian can be diagonalized ex- 
actly for small systems, using the Lanczos method. 47 ' 48 
Unless otherwise noted, the results presented here are 
on a 4 x 3 unit-cell torus (i.e. 4x6 real-space sites). 
We do not project onto the flat subband, see Fig. 4(a), 
but instead model the whole system of both subbands 
with C = ±1. While this approach increases the Hilbcrt 
space and thus restricts system size, it has the advan- 
tage that competition with phases mixing the two sub- 
bands is included automatically. We use several observ- 
ables to detect FCI states and to distinguish them from 
other phases: first, we discuss conclusions to be drawn 
from the eigenvalues, and later also include information 
obtained from the eigenstates, namely the many-body 
Chcrn number as a topological invariant and the charge- 
structure factor indicating formation of a (conventional) 
charge-density wave. 

Even though the results of our investigation do not 
allow for conclusive finite-size scaling, we have verified 
that the findings presented here are consistent for smaller 
as well as somewhat larger systems. Furthermore, the 
main features of the FCI states found are consistent for 
all aspect ratios yielding the same system size, apart from 
those that reduce the system to a one-dimensional chain. 
Before passing, we note that we have obtained similar 
results to the ones presented here for several other filling 
fractions v = p/q, with q = 5,7. The states at filling 
fractions v = 4/5,5/7,6/7 are however elusive, as has 
also been found to be the case in another model. 24 



A. Eigenvalues and flux insertion 

FCI states manifest themselves in features of the ob- 
tained eigenvalue spectra. Traditional FQH ground 
states are g-fold quasi-degenerate on a finite torus at fill- 
ing v = p/q. i9 The same degeneracy is expected to occur 
for FCI states as well, at least within certain well-defined 
limits. 30 Quasi-degenerate FQH ground-state eigenvalues 
must also exhibit spectral flow, leading each of them 
into another upon insertion of a flux quantum. 50 Mag- 
netic fluxes can be modelled by introducing phase fac- 
tors to the hopping from site i = i x a\ + iy&2 to site 
j = j x ai + j y SL2, thus leading to the transformation 
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FIG. 5. (Color online) Eigenvalue spectrum and evolution of 
selected levels under flux insertion for different values of V/t. 
In all panels t' /t = 0.2. 



tion picks up a phase 4> x / y . 

The features described above are illustrated, as a func- 
tion of V/t, in Fig. 5 for v = 1/3. In this case, the three 
quasi-degenerate ground states emerge from the contin- 
uum as the interaction strength is increased. Insertion 
of one flux quantum indeed leads from one ground-state 
eigenvalue to another. It can also be seen that spectral 
flow is a general property of the eigenvalue spectrum, 
and we have found that it also occurs even at V/t = 0. 
The v = 2/3 ground states behave in a similar fashion. 
It can be seen that the FCI states remain gapped for a 
range of both M and V. In order to define a physically 
meaningful phase boundary however, one has to make 
sure that the ground states remain gapped upon flux in- 
sertion, not only in their own momentum sector, but also 
across momentum sectors. The reason is that the pres- 
ence of impurities would mix the momentum sectors and 
therefore close the gap if the quasi-degenerate FCI-state 
levels cross excited-state levels. 

As seen in Fig. 5, interaction mostly changes the en- 
ergy differences between groups of three quasi-degenerate 
eigenstates, which flow into each other upon flux inser- 
tion, and has less impact on the energy splitting within 
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FIG. 6. (Color online) Eigenvalue spectrum and evolution of 
selected levels under flux insertion for different values of t'/t. 
In all panels V/t = 1. 



each manifold. When tuning the transition via increasing 
the band dispersion, the situation is somewhat different, 
see Fig. 6. For M « 10, the FCI ground state is de- 
stroyed by increasing the split between quasi-degenerate 
levels, permitting a higher-energy states to mix with 
the FCI-state manifold for some flux values. The split- 
ting between ground states in the traditional FQHE case 
is due to quasiparticlc-quasiholc excitations propagating 
around the torus and therefore depends mainly on sys- 
tem size, becoming zero in the thermodynamic limit. 49 
In the finite tight-binding lattice systems discussed here, 
the splitting is affected significantly by the residual ki- 
netic energy of the partially filled band, as can be seen 
in Fig. 6. Since the dispersion is present also in the ther- 
modynamic limit, it is not clear whether the splitting 
survives to the thermodynamic limit or not in this sys- 
tem. 

Despite level rearrangements or increasing spread be- 
tween eigenvalues, FCI states remain topologically con- 
jugate upon varying the band flatness or the interaction 
strength. A factor that can split this conjugacy is dis- 
order, as will be shown in the Sec. HID. At the same 
time, ground-state levels that remain gapped upon flux 
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FIG. 7. (Color online) Eigenenergy spectrum and spectral 
flow for FCI state with quasiholes introduced (a) by enlarging 
the system or (b) by removing particles. The numbers of 
levels below the legends are (a) 35 and (b) 40 respectively, in 
agreement with the counting rule of Ref. 30. In the 7x2- 
unit cell system, the filling fraction is v — 2/7. There are 
seven quasi-degenerate ground states (marked in red), which 
are slightly separated from the rest of the states in the low- 
energy sector and exhibit spectral flow (shown in different 
colors). Their Chern number is C = 2/7 within the numerical 
error margin. 



insertion along one direction may cross excited-state lev- 
els upon flux insertion along another direction. A more 
careful analysis, based both on eigenenergy as well as on 
cigenstate properties, is necessary to uniquely determine 
a gapped ground state for a given parameter set. This is 
the purpose of Sec. Ill C. Before investigating eigenstates, 
however, the next section is going to discuss another as- 
pect of eigenenergies, fractional statistics. 



B. Fractional statistics 

The concept of generalized Pauli principles, 6 accord- 
ing to which states with clustered anyonic quasiparticlcs 
are energetically penalized, has served as a tool to indi- 
cate fractional statistics in FQH states. 51,52 As was first 
demonstrated heuristically 20 and later on substantiated 
theoretically, 30 the same logic holds for FCI states. For 
the v = 1/3 FCI state, the number of quasihole-state 
levels in a well separated, low-energy Fock space is equal 
to the number of (l,3)-admissiblc partitions of the mo- 
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mentum sectors on the torus. The individual states per 
momentum sector can also be accounted for, due to an 
emergent translational symmetry (sec Rcf . 30 for details) . 
This symmetry is a characteristic feature of all the FCI 
states we have found in the triangular lattice model. In 
this section, wc demonstrate that quasihole states of the 
triangular-lattice model obey the same state-counting ar- 
guments. 

Quasiholes can be introduced in a FCI state by either 
removing electrons or increasing the system size. An ex- 
ample of both cases is presented in Fig. 7. The count- 
ing rule is verified in both cases. The counting rule can 
also be easily generalized to all Laughlin filling fractions 
v = 1/q and wc have verified it for v = 1/5 in this model. 

Apart from the agreement to the counting rule, one 
can also notice the emergence of a "daughter" FCI state 
within the gapped low-energy sector, as expected in the 
hierarchy picture. Even though in the present context 
this is mainly a peculiarity of the small system sizes, it 
can nevertheless be viewed as a simple example of the 
formation of a hierarchical FCI state among the quasi- 
particle states of a "parent" FCI state. Taking this ob- 
servation one step further, it can be seen that the spec- 
trum looks qualitatively different under flux insertion in 
the two cases presented in Fig. 7, namely, the eigenvalue 
spectrum in the case where the filling is v = 1/4 looks 
more like that of the non-interacting system with levels 
crossing upon flux insertion, whereas at v = 2/7 lev- 
els come in small groups, which remain separated un- 
der flux insertion (compare to Fig. 5). This observation 
can be compared to the composite fermion theory of the 
FQHE, 53 according to which the state of a FQH sys- 
tem at filling fractions with even denominators can be 
effectively described by free fermions with flux tubes at- 
tached to them. This compositc-fermion view of the FCI 
hierarchy of states has also been supported by recent nu- 
merical calculations. 24 These two different facets of the 
same model can serve as an example of the compatibility 
between the hierarchy and composite-fermion pictures of 
the FQHE in an unconventional FQH system. 




CM 

ft 



12 3' 100 200 300 




FIG. 8. (Color online) (a-c) Berry curvatures, (d-f) their stan- 
dard deviations from the exact average and relative deviations 
of Chern numbers as a function of grid size for the same state, 
when it is in the excited-state quasi-continuum (top) and in 
the ground-state manifold (middle and bottom). The integra- 
tions in Eq. (18) have been approximated by simple Riemann 
sums. N$ is the number of <j> x / y points taken in the range 
[0, 6n) in each direction. 



C. Topological invariant 

The most unambiguous characteristic of a FCI state 
is arguably its Chern number. The ground-state Chern 
number of a FCI at filling v = p/q should be exactly v. 
The Chern number of a gapped, degenerate state can be 
efficiently calculated using the Kubo formula: 54 ' 55 

(18) 

where N c is the number of unit cells, |n) are the g-fold de- 
generate ground states and \n'} are higher-energy eigen- 
states. The corresponding eigenenergies are e n / n i. The 
^/^-dcpcndcnce of the Hamiltonian comes from mag- 



netic fluxes going through each of the handles of the 
torus, as discussed in Sec. II C. 

The integrand in Eq. 18 is proportional to the Berry 
curvature B n for each of the states in the degenerate 
ground state. For FCI states, Berry curvatures are pe- 
riodic functions varying with <f> x /y When the single- 
particle Berry curvature of a Chern band is smooth, cer- 
tain theoretical tools can be exploited in order to draw 
analogies with the corresponding FQHE state. 30,37,56 In 
Fig. 8, the many-body Berry curvature for a specific 
state, which develops into one of the FCI ground states, is 
shown. It is seen that its smoothness varies as this state 
emerges from the excited state continuum. The period 
extends over q flux quanta in one direction and one flux 
quantum in the other, remains unchanged for all values 
of V and it is the same for all states in the ground-state 
manifold. The Berry curvature of the other two ground 
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FIG. 9. (Color online) Chern number C, see Eq. (18), for 
the three quasi-degenerate FCI states as a function of (a) the 
interaction strength V/t and (b) the third-neighbor hopping 
t'/t (right). Not the different scales on the vertical axes. 



states is the same function, but translated by one and 
two flux quanta in the <p v direction respectively. Even 
when the many-body Berry curvature of a FCI state is 
strongly varying, it is often centered around the value 
corresponding to the filling fraction and, since it is peri- 
odic, the Chern number obtained by integration is close 
to the quantized value. This occurs even if the FCI state 
is not a ground state. In the example presented in Fig. 8, 
the Chern number is in all three cases equal to v = 1/3 
within the numerical accuracy, which we will now discuss. 

The Chern numbers we have calculated are very close 
to the expected exact values in the FCI regime, despite 
the fact that the integrations involved are performed nu- 
merically. The relative deviation with respect to the ex- 
act value for the case of v = 1/3 is presented in Fig. 8. 
This example illustrates that the error in the Chern num- 
ber due to the finite size of the system should be smaller 
than 1%, at least for the Kubo formula approach, even for 
small systems. In the case of clearly gapped FCI ground 
states, the accuracy is even better, due to the integrand 
being very smooth. In all cases, we have used simple Rie- 
mann summation to evaluate integrals, in order to obtain 
upper bounds for the numerical errors. Other methods, 
like the Simpson rule, would converge for smaller grid 
sizes. It should be noted that when the Chern num- 
ber of an individual state within the degenerate ground 
state is evaluated, the integration should be extended 
over the whole period of the Berry curvature. This is 
no longer necessary when calculating the average Chern 
number and the integration range can then be restricted 
to [0,2tt) 2 . 

The effects of interaction and band dispersion can also 
be traced in the behavior of the Chern number. This 
is shown in Fig. 9. Outside the FCI regime, where low- 
energy levels cross upon flux insertion, the Chern number 
oscillates. It converges to v at about the value of V/t 
for which the FCI states separate from the excited-state 
spectrum, forming a three-fold degenerate ground state. 
On the other hand, the Chern number of these states 
remains close to 1/3 with deviation from the fiat-band 
limit, and only changes smoothly with t'/t, converging to 
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FIG. 10. (Color online) (a,b) Chern numbers C as a function 
of the disorder phase <p. (c,d) Standard deviation of total 
Berry curvature B, defined by adding up the many-body Berry 
curvatures B n of the degenerate ground states for each value 
of 4>x/ y , and relative deviation of Chern number as a function 



of grid size. N<j, is the number of 



c/a 



points taken in the 



range [0, 2n) in each direction. The results shown in the right 
panels are for <p ~ 2n/3, but are qualitatively the same for all 
other values of (p. 



the expected value close to the phase boundary discussed 
in Sec. Ill E. 



D. Disorder 

FQH states must be robust against disorder, as long 
as the energy scale of the disorder is smaller than the 
gap. 49 The magnetic texture in the three-orbital model 
presented is generated spontaneously by itinerant elec- 
trons, so an expected source of disorder are local inho- 
mogencities in the chiral spin pattern. A simplistic ap- 
proach to simulate this effect is to vary the flux picked 
up by an electron hopping along one selected bond of the 
finite cluster. This is done by varying the phase tp in the 
phase factor in front of the corresponding hopping. The 
effect of such a variation on Chern numbers at v = 1/3 
is shown in Fig. 10. 

As long as the energy scale of the disorder is small 
enough, the ground states remain separated by a gap 
from the rest of the spectrum and their average Chern 
number remains constant. The Chern numbers of the in- 
dividual quasi-degenerate ground states are not smooth 
functions of magnetic disorder. In particular, as soon as 
the impurity is switched on, the Chern number of each 
individual state jumps to an integer value. Also, Berry 
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curvatures are no longer smooth functions of 4> x / y . Nev- 
ertheless, the Chern numbers are still numerically well- 
defined, as demonstrated in the right panels of Fig. 10, 
and their average, which is the proper observable quan- 
tity in the thermodynamic limit, remains constant as the 
disorder is varied. This invariance of the Chern number 
directly demonstrates the topological robustness of FCI 
states. 

Again in analogy to FQH states, 57 disorder introduces 
a splitting between FCI ground-state eigenvalues, which 
no longer exchange places under flux insertion. This 
splitting is expected to be a finite-size effect and should 
disappear in the thermodynamic limit. A conclusive 
proof of this statement is however beyond the scope of 
this work. Despite the similarities, the impact of disor- 
der in general on the energy-scale balance necessary for 
FCI states is qualitatively different than the situation in 
traditional FQH states. For example, disorder may have 
a different effect depending on the interaction strength 
and range, and may lead to another, possibly topologi- 
cally trivial, state, as has been shown in a recent study 
on the checkerboard lattice involving chemical-potential 
disorder. 58 



E. Competition with Charge-density wave and 
Phase diagram 




FIG. 11. (color online) Dispersion (gray-/ color-scale) and 
Fermi surface (FS) for v = 2n = 2/3 (thick solid line) for (a) 
t'/t = 0.16 and M = 11.16, (b) t'/t = 0.19245 and M = 24 
(very close to maximal M) and (c) t'/t = 0.25 and M = 5.5. 
Dashed lines indicate the first Brillouin zone, (d) illustrates 
the CDW possible for v = In = 2/3. Black and red/gray 
circles indicate the two sublattices of the spin-chiral order 
and the effective model Eq. (15), filled circles the particles in 
the CDW state induced by large V/t. 



Early on, it was pointed out that an important prereq- 
uisite for the emergence of a FCI state is the balance be- 
tween three energy scales, namely the width of the topo- 
logically non-trivial band, the gap(s) separating it from 
other bands and the Coulomb interaction strength. If the 
ratio between the first two, as expressed by the figure of 
merit Eq. (5), is large, then the interaction can become 
strong enough to induce FCI states, but remain small 
compared to the band gaps, which avoids mixing in wave 
function of different topological character. Many inves- 
tigations have been focused on models in the perfectly 
flat band limit (exactly as in a Landau level), supple- 
mented by various types of interaction. Rather recently, 
it was reported that the band flatness alone is not in fact 
a reliable indicator for the stability of FCI states. 29 

In this section, we discuss how a finite dispersion in- 
fluences the stability of FCI states in a triangular-lattice 
model, where we focus on band filling and FS nesting. 
The first presents a rather obvious difference to LLs, as 
has recently also been mentioned in Ref. 24: While a 
LL is expected to be particle-hole symmetric, the nearly- 
flat subband of a lattice model is not. We are going to 
discuss filling fractions v = 1/3 and v = 2/3 that turn 
out to clearly exemplify this difference. The latter case, 
which was also shortly discussed in Ref. 22, corresponds 
to 1/3 filling of the original triangular lattice, where NN 
Coulomb interaction can stabilize a CDW [with particles 
at second-neighbor sites, see Fig. 11(d)], while no such 
CDW is possible at v = 1/3. The competition between 
the FCI and the CDW is in turn strongly affected by FS 



nesting, and we are going to see that the FCI is far more 
stable for more dispersive bands without nesting than for 
flatter bands with a nested Fermi surface. However, we 
are also going to see that very flat bands allow for FCI 
states at the lowest interaction strengths, even though 
they are in our case very well nested. 

As mentioned in Sec. II C, varying t'/t allows us to 
tune the band flatness and to moreover switch between 
regimes with and without FS nesting. Examples are 
shown in Fig. 11, where the dispersion as well as the 
Fermi surface corresponding to v = 2/3 are shown for 
some values of t'/t. Both for near-optimal t'/t = 0.19245 
(M = 24) and for smaller t'/t = 0.16 (M « 11), the 
FS contains hexagons with almost perfectly nested seg- 
ments. A difference between the two cases is on one hand 
the flatness ratio, but on the other hand, the flatter bands 
also have an additional circular FS around the T point. 
For t'/t > 0.23, only the circular FS remains (see the 
example with t'/t = 0.25 and M = 5.5) and there is thus 
no longer good nesting. 

Figure 12 shows the region where FCI states are sta- 
ble on a 24-site (4x3 unit-cell) system, as a function 
of Coulomb interaction strength V and the flatness ratio 
M, which is in turn controlled by third- neighbor hop- 
ping, see Fig. 4. The two panels are for fillings v = 1/3 
and v = 2/3. In both cases, the ground-state manifold 
in the FCI state is expected to be three-fold degenerate. 
The colorcode indicates the gap between the three ground 
states and the fourth lowest eigenstate in the absence of 
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FIG. 12. (Color online) Gap of three-fold degenerate FCI 
ground state as a function of M and V/t, for (a) v = 1/3 and 
(b) v — 2/3, shown in colorcode. The dashed lines indicate 
the phase boundary, defined by the condition that the gap 
remains open for all values of inserted magnetic flux. The 
flatness ratio M of the flat band of //kin (bottom scale) is 
adjusted by varying t'/t (top scale). The maximum value of 
the flatness ratio (M ~ 24) is marked by the dotted lines. 



applied magnetic flux. In order to define a physically 
meaningful phase boundary however, one has to make 
sure that the ground states remain gapped upon flux in- 
sertion, not only in their own momentum sector, but also 
across momentum sectors, see Sec. Ill A. The solid lines 
in Fig. 12 indicate the phase boundary determined by 
taking these considerations into account. We have veri- 
fied that the phase diagrams obtained from a 30-site torus 
arc qualitatively the same as the ones shown here. The 
stability range of the FCI differs for v — 2/3 and system 
sizes commensurate with the CDW (3x6 and 6x6 sites) , 
but the corresponding phase-diagram area is still finite. 

For both fillings, the system is in a metallic state at 
small V/t, while finite values of V/t can lead to FCI 
ground states for an extended range of the band flatness. 



FIG. 13. (Color online) Static structure factor for various 
values of t'/t and V/t at filling fractions (a,b) v = 1/3 and 
(c-f) v = 2/3. The sharp peaks indicate the formation of a 
charge-ordered state. 



In the case of v = 1/3, the FCI state can be induced for 
any band dispersion by making V large enough. The FCI 
persists even for V/t considerably larger than the band 
gap. This has also been observed for another FCI model 
on the checkerboard lattice, 18 and appears to indicate 
that the fermions are dilute enough to occupy mostly 
one subband regardless of V/t, so that a mixing of the 
topological character of the two bands does not occur. In 
contrast, the v = 2/3 FCI states survive only for moder- 
ate values of V/t. Larger interaction strengths lead to a 
different ground state, which is, as we will argue below, 
a CDW. 

To determine the type of order in each of the regions in 
Fig. 12, we calculate key properties of the corresponding 
ground states. The FCI states, despite the hints from the 
energy spectrum, have to be identified by their topolog- 
ical invariant, see Sec. IIIC. To find out whether there 
is tendency towards charge order in the rest of the phase 
diagram, we calculate the static (charge-)structure factor 
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(SSF), defined as: 
1 Ns 

»(k) = ^(0|E e ik - (R ^ Ri) (^ - »)(n, - n)|0) , (19) 

3,1 

where N s is the number of sites, fij is the electron-number 
operator acting on the site at position Rj, n = v/2 is the 
average electron number and |0) stands for the ground 
state. Charge-density modulations are marked by sharp 
features in n(k) at certain wavevectors. Liquid states, 
on the other hand, should be featureless in comparison 
to charge-modulated states. 

Even though the accessible system sizes are not 
large enough to exemplify the featurelessness of liquid 
states, a qualitative difference between liquid and charge- 
modulated states can be seen. Examples are shown in 
Fig. 13. At v = 1/3 and within the FCI regime, n(k) 
remains almost unchanged upon variation of model pa- 
rameters. The same holds, although less markedly, for 
the metallic, Fermi-liquid-like state at small interaction 
strengths and small (or very large) t'/t. The shape of 
n(k) is nevertheless distinct for the two different liq- 
uid states, see Figs. ll(d)(a) and (b). At v = 2/3 and 
t'/t = 0.1 < 0.18, away from the FCI regime, two peaks 
appear in n(k) and increase continuously upon increasing 
V/t, see see Fig. 13(d). When choosing t'/t = 0.2 near 
the maximal flatness [Fig. 13(c)] or t'/t = 0.35 with bad 
FS nesting [Fig. 13(f)], the peaks only begin to grow for 
large V , where ground-state behavior upon flux insertion 
changes: a full gap is only obtained for lattice sizes com- 
mensurate with the ordering pattern, like the 3x6 and 
6x6 lattices, but not on more general lattices. This indi- 
cates that the FCI state breaks down and is replaced by a 
CDW. The peaks observed in n(k) for v = 2/3 grow with 
V, supporting their relation to a CDW. Their wavevec- 
tors correspond to a state where particles sit at second 
nearest-neighbor sites on the triangular lattice, i.e., to 
the regular charge pattern compatible with a filling of 
n = v/2 = 1/3 of the triangular lattice. As can be seen 
in Fig. 11(d), this pattern avoids any penalty due to NN 
Coulomb interaction V and large enough V ^> t, t' will 
thus eventually induce such a charge distribution. 

The phase diagrams in Fig. 12 clearly show that the 
band flatness parametrized by M is itself not a reliable 
indicator for the stability of FCI states, as has also been 
pointed out recently in a different case. 29 The energy gap 
for v = 1/3 is rather symmetric with respect to the high- 
est figure of merit, but the phase boundary determined 
by requiring a full gap over all momentum sectors and for 
all fluxes shows that the FCI states are somewhat more 
stable for larger t'/t. In the case of v = 2/3, the asym- 
metry is far more striking: as one can see in Fig. 12(b), 
FCI states require rather flat bands with M 13 for 
t'/t < 0.2, but extend to a band with a width compa- 
rable to the gap separating it from its counterpart for 
t'/t > 0.2. Having rather flat bands indeed makes it eas- 
ier for small V/t to induce FCI states both at v = 1/3 
and v = 2/3, even though the optimal M is still not quite 



the largest, at least for our system sizes. As soon as the 
bands acquire some dispersion, however, features beyond 
band flatness, in our case FS nesting, can strongly in- 
fluence the stability of FCI states by favoring competing 
states, in our case a CDW. Another key feature, which 
extends previous results, 29 is that in the case of competi- 
tion between FCI and other phases, perfect band flatness 
is not necessarily the ideal condition for the stability of 
FCI ground states. 

Before passing, a few more comments on the phase 
diagram have to be made. Ideally, a ground-state prop- 
erty would be used to determine the phase boundaries. 
One such property, which is sensitive to phase transi- 
tions, is the ground-state fidelity, defined as a measure 
of the overlap (Tp(a)\ip(a + 6a)}, where ip is the ground- 
state wavefunction and a is a control parameter varied in 
small steps 5a. Phase transitions are then marked by a 
divergence in the fidelity at the transition point. Having 
calculated the ground-state fidelity upon varying inter- 
action strength and band flatness for the cases presented 
here, we find that such divergences occur only at points 
where ground-state level crossings also occur, so the fi- 
delity does not provide any extra information compared 
to the eigenvalue spectra. Furthermore, it has recently 
been shown that the fidelity on finite systems can fail 
to register topological phase transitions, e.g., between a 
FQH-like state and a Fermi liquid. 58 



IV. SUMMARY AND CONCLUSIONS 

Based on an earlier investigation, 22 we have illustrated 
extensively how FCI states can emerge in a strongly cor- 
related multiorbital model and have shown how an ef- 
fective spinless one-orbital model with nearly flat bands 
and non-zero Chern number C = 1 arises as the the low- 
energy limit of Kondo-lattice and Hubbard models for 
t2 g orbitals on a triangular lattice. NN Coulomb interac- 
tion V then stabilizes states with all the characteristics of 
FCI states: The lowest eigenvalues have a near degener- 
acy corresponding to the denominator of the filling frac- 
tion and show spectral flow without closing of the energy 
gap. We have moreover demonstrated that the states of 
the ground-state manifold of this phase have a non-zero 
Chern number, which is shown to be precisely quantized 
and equal to the filling fraction. The exact quantization 
of the Chern number holds also when FCI states are not 
ground states, and survives when the many-body Berry 
curvature is not a particularly flat function. 

Deviating from the exact filling fraction v = 1/3 by ei- 
ther removing electrons or increasing the system size, we 
found indications of fractional quasiholc statistics in the 
eigenvalue spectra, by application of the state-counting 
rule elaborated in Ref. 30. The eigenvalue spectra con- 
tain features which remind one of both the hierarchy and 
the composite fermion pictures of FQH states, therefore 
pushing the correspondence between FCI and FQH states 
one step further. 
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Since in the present context the magnetic texture is 
generated by strongly correlated itinerant electrons, im- 
perfections in this texture are to be expected. We there- 
fore investigate the impact of a localized impurity that 
mixes states across momentum sectors. Disorder causes 
splitting of the energy levels corresponding to the degen- 
erate FCI ground states, thus making these states in- 
equivalent in finite systems. However, the gap between 
these states and the excited-state spectrum remains open 
and their average Chern number remains accurately fixed 
to the value of the filling fraction, even though the Chern 
numbers of the individual states deviate from that value. 
This behavior illustrates the topological protection of the 
FCI states. 

Finally, we address the robustness of FCI states with 
respect to band flatness, interaction strength, filling and 
FS nesting. Very flat bands indeed turn out to require 
the weakest interaction to produce an FCI state, corrob- 
orating the figure of merit M [see Eq. (5)] as a good first 
step in assessing candidate systems. As soon as one de- 
viates from this limit, however, additional properties of 
the system become crucial: In contrast to a perfectly flat 



LL, strong particle-hole asymmetry can arise on a lat- 
tice, because some fillings allow CDWs more easily than 
others. In addition, a CDW is favored by FS nesting 
and this in turn reduces the stability range of the FCI 
states. Nevertheless, it turns out that FCI states can be 
obtained if either of the following conditions is met: (i) 
very fiat bands, even near nesting and at a filling favor- 
able to a CDW; (ii) even rather dispersive bands, if the 
filling fraction or the absence of FS nesting are unfavor- 
able to a CDW. This flexibility makes FCI states appear 
more realistic. 
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